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The Hartree-Fock approximation to the many-fermion problem can break exact symmetries, and in 
some cases by changing a parameter in the interaction one can drive the Hartree-Fock minimum from 
a symmetry-breaking state to a symmetry-conserving state (also referred to as a "phase transition" 
in the literature). The order of the transition is important when one applies the random phase 
approximation (RPA) to the of the Hartree-Fock wavefunction: if first order, RPA is stable through 
the transition, but if second-order, then the RPA amplitudes become large and lead to unphysical 
results. The latter is known as "collapse" of the RPA. While the difference between first- and 
second-order transitions in the RPA was first pointed out by Thouless, we present for the first time 
non-trivial examples of both first- and second-order transitions in a uniform model, the interacting 
shell-model, where we can compare to exact numerical results. 

I. INTRODUCTION 

An important class of approximations in many-body theory are mean- field approximations, which reduce the many- 
body problem to an effective one-body problem, and their generalizations [ij. In this paper we consider the Hartree- 
Fock approximation, a variational approach that approximates the ground state wavefunction by a single Slater de- 
terminant (antisymmeterized product of single-particle wavefunctions), and the random phase approximation (RPA), 
which builds small amplitude correlations on top of the Hartree-Fock state, and which itself can be derived as the 
small-amplitude limit of the time- dependent Hartree-Fock approximation. 

Because the Hartree-Fock (HF) state ignores correlations, it can break exact symmetries, such as translational 
and rotational invariance. A related approximation, the Hartree-Fock-Bogoliubov (HFB) approximation, built upon 
quasi-particles, also breaks conservation of the number of particles. Despite this breaking of symmetries, HF, HFB, 
' ■ 1 ' and the RPA and quasi-particle RPA are widely and successfully used to describe low-energy nuclear structure. 

, If one dials the parameters of the Hamiltonian, it is possible for the ground state HF or HFB wavefunction to 
change from a symmetry-preserving state to one that breaks an exact symmetry. For example, as one changes the 
single-particle energies as one attempts to fit data, the HF state may go between spherical and deformed states, and 
the HFB state may between normal (number-conserving) and supcrfluid (number-nonconserving) states. 

A long time ago Thouless pointed out 0] that, as the mean- field state is driven between symmetry-conserving 
(SC) and symmetry-nonconserving (SNC) states, there are two possible kinds of state transitions (often referred to 
as "phase transitions" in the literature, even for finite systems). For second-order transitions the fluctuations in RPA 
become unphysically large, leading to the so-called "collapse" of the RPA. A frequently cited example is the collapse 
5_i ■ of RPA in the Lipkin-Meshkov-Glick model Significant effort in the literature has been devoted to collapse in the 
RPA and QRPA and possible solutions. All but forgotten are first-order transitions for which RPA does not collapse. 

In this paper we use a non-trivial realization of the HF and RPA in the interacting shell model, where we can 
compare to exact numerical calculations, and demonstrate both first- and second-order transitions. As predicted 
by Thouless, second-order transitions and collapse are driven by odd-parity operators or modes, while first-order 
transitions are driven by even-parity operators or modes. In this paper we do not propose any new solutions to the 
collapse of RPA at or near second-order transitions; instead, this work illustrates that the world of RPA and state 
transitions is more complicated than the standard narrative. 
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II. HARTREE-FOCK AND RPA CALCULATIONS IN THE INTERACTING SHELL MODEL 

We work in the framework of the interacting shell model. The model space is defined by a truncated set of single- 
particle orbits, for example the 0pi/ 2 -0p3/2 space, usually called the p-shell, or the lsi/2-0c?3/2-0(i5/2 space, the sd 
shell. The interaction is given by a set of single-particle energies e a and two-body matrix elements VjT(ab,cd). The 
interactions used are not necessarily simple schematic forces but generally start from carefully computed G-matrix 
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interactions and then adjusted empirically to reproduce a large number of ground state binding energies and excitation 
energies. 

There are a number of computer programs that read in the model space, generate many-body basis states (Slater 
determinants in occupation space), compute the many-body matrix elements from the single-particle energies and 
two-body matrix elements, and diagonalize the Hamiltonian matrix to get out the eigenenergies and wavefunctions, 
from which one can compute observables, transitions, etc.. For this paper we used the REDSTICK shell model code||. 

For Hartree-Fock plus RPA, we have use the SHERPA (SHEll-model RPA) code [6[ which uses exactly the same 
input as REDSTICK. SHERPA places no restriction on the Hartree-Fock wavefunction except that it must be purely 
real; otherwise it can have any arbitrary deformation contained in the model space. Thus one can have spherical, 
axially symmetric deformation, triaxial deformation, and even parity-mixed ground state. 

In previous papers we have used SHERPA to directly test HF+RPA as an approximation to full shell-model 
diagonalization, looking at correlation energies B, ground state observables Q, electromagnetic transitions 0, and 
charge-changing Gamow- Teller transitions [10(. In presenting those results we were frequenly asked about the issue of 
the collapse of RPA. Such questions motivated this paper. 



III. FIRST- AND SECOND-ORDER TRANSITIONS 

In this section we revisit Thouless' arguments 0, Q regarding the order of the transition. To do this, we need to 
remind the reader of at least one way to develop RPA fij , although we do not give the derivation in full. 

Let |^/} be a Slater determinant, that is, an antisymmeterized product of single-particle wavefunctions. In second 
quantization, where a\ creates a particle in state a, we write 

i*)=n^i°)- 

i 

We follow the usual convention where i,j denote occupied states and m, n denote unoccuped states. One then 
introduces the general particle- hole operator, Z^ = ^2, ni z^a^a^, (where the z m i can be complex) and the state 

\z) =exp(Z)|*> (2) 
which, by Thouless' theorem Q is a Slater determinant not orthogonal to the starting state. One can compute 

'zH z" 



E{z) 



(z\z) 



(3) 



which, assuming z small, can be expanded 



E(z) — Eq + ^ ' h mi Z m i + h m iZ mi + ^ ' ^-mi,nj z mi z n j + ^ Bmi,nj :Z mi Z nj 2 ^mi,nj z mi z nj + (4) 

mi mi.nj 

At a local minimum, h m i — 0; this is the Hartree-Fock condition, and Eq is the Hartree-Fock energy. The quadratic 
terms can be treated as a harmonic oscillator: one treats the z mi as boson operators and using a Bogoliubov trans- 
formation put into diagonal form, using the famous RPA matrix equation: 

A B \( X x \ _ Q f X X \ (5) 



B* -A* J \Y X J V Y x 

There can be, of course, higher-order corrections in the energy landscape beyond quadratic, which we will consider 
shortly. 

Suppose one is in a SC state, e.g., a state of good angular momentum, usually spherical, or parity. Then, solving 
Eq. ([5]) one finds the X, Y modes also have good symmetry. For example, with our code SHERPA, if the HF state is 
spherically symmetric, one gets out RPA modes that have good angular momentum J, and one sees the appropriate 
(2J + 1) degeneracy in the RPA frequencies (If, on the other hand, rotational symmetry is broken but one has 
axial symmetry, then one can have two-fold degeneracies, signaling RPA modes that are time-reversed of each other, 
or a single mode, which must be time-reversal even. In the event of triaxiality, one has no degeneracies in the RPA 
spectrum.) 

To consider "state transitions," from SC to SNC, one needs to look at higher-order terms. Let Z\ represent 
generically the RPA modes; by expanding about the HF minimum, one expands the energy landscape by (HF\Z™\HF). 



FIG. 1: Sketches of "phase" transitions in the energy landscape, (a) A first-order transition, (b) A second-order transition. 



Because we are at a minimum, n = 1 vanishes, n = 2 yield the curvature and the RPA frequencies. What about 
71 = 3,4,...? 

Now we get to the heart of Thouless' argument. If the RPA mode Z\ has odd parity, while the HF state has good 
parity, then (HF\Z™\HF) must vanish for all odd n. On the other hand, if Z\ has even parity, then (HF\Z\\HF) 
can be nonzero; for example, it is possible to couple three J ~ 2 operators to total angular momentum zero. 

Figure Q] illustrates both cases. Fig. QJa) is f° r modes with even parity, so that cubic terms play a role. One can 
clearly have coexisting local minima, and so one gets a first-order transition. Fig. [Hb) is for odd-parity modes, so 
that the energy landscape must be symmetric. One tends to get only a second-order transition. 

In our examples below, we see Thouless' predictions played out. A system with only even parity modes has first-order 
transitions, while a system with odd parity modes has a second-order transition. 

What happens in a second-order transition? In that case the RPA frequencies fl\ — > and the fluctuations about the 
HF state, measured by \Y\\, become very large, allowed because the RPA vectors have a nonstandard normalization: 
|A"a| 2 - \Y\\ 2 = 1. The RPA energy is 

E RPA = E RF - Y,M\\Yx\ 2 ~ tItt ( 6 ) 
I^o 2M ° 

(see 0, 0] for details), and so as |Y\| 2 becomes large, the energy dives or collapses. The appearance of unphysical 
values is unsurprising because in the derivation one assumes small amplitudes for X and Y. 

In the published literature, discussions regarding the behavior of the RPA near a phase transition focus exclusively 
on collapse of RPA 0, 

Hlj |. that is, on second-order transitions. Outside of Thouless 0] there is no discussion of 
first-order transitions (indeed, the literature appears to uniformly refer to phase transitions without quantifying the 
order of the transition). Yet, as illustrated below, it is not hard to find a first-order transition. It is possible that 
people using RPA have encountered first-order transitions without realizing it, as there is no catastrophic collapse to 
signal the transition. While clearly second-order transitions are more problematic, we find it instructive to explore 
both kinds. 



IV. RESULTS 

A. Example of a first-order transition: deformation 

We begin with two case studies in the sd-shell, using Wildenthal's universal sd (USD) interaction[12]. The first of 
these is 28 Si (6 valence protons and 6 valence neutrons). By increasing the difference A between the 0d 5 / 2 single- 
particle energy and the lsi/ 2 -0d 3 /2 single-particle energies, we can force the Hartree-Fock state to go from an oblate 
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FIG. 2: (Color online) Ground state energy of 28 Si, calculated in the full interacting shell model (SM, solid lines), Hartree-Fock 
(HF, dotted lines), and Hartree-Fock plus random phase approximation (HF+RPA, dashed lines). Here A is added to the lsi/2 
and Qd 3 /2 single-particle energies; A = corresponds to the original Wildenthal values. The left-hand panel shows only the 
final result; the right-hand panel identifies spherical and (oblate) deformed mean- field phases. Note the significant region of 
coexistence. 



deformed state to a spherical state with the shell filled. For convenience A = corresponds to the original 

Wildenthal values. 

Fig. @ shows the "exact" calculation, which is an interacting shell-model calculation performed in the full Ofiuj 
sd valence space, compared to the lowest Hartree-Fock energy, and the RPA correlation correction on top of the HF 
energy. Because one switches between two degenerate HF states, the HF energy is continuous, while the HF+RPA 
energy is discontinuous, because the curvatures (RPA frequencies) are different. The right-hand panel illuminates in 
detail what is going on. Here we explicitly show the (oblate) deformed and spherical HF energies and their respective 
RPA corrections. One sees there is a significant region of parameter space, encompassing the original Wildenthal 
value, where locally stable deformed and spherical HF solution coexist. 

As we drive A further positive or negative, eventually the deformed or spherical solutions, respectively, become 
unstable. In Fig. ^ this is seen as the HF+RPA energy dives sharply. For further diagnosis, in the upper half of 
Fig. ([3]) we plot the lowest nonzero RPA eigcnfrcqucncics £1 as a function of A. (The deformed state also has two 
zero- frequency modes corresponding to broken rotational symmetries.) We see that one eigenfrequency dives to zero, 
signalling an instability. As further diagnosis, we give the degeneracy of the RPA eigenfrequencies; for the spherical HF 
state, the degeneracy of the collapsing eigenfrequency is 5, suggesting a quadrupole mode. For the oblate deformed 
state, the RPA eigenmodes either come in time-reversed pairs (degeneracy = 2) or are already time-reversal-even 
(degeneracy =1). 

Note: We cannot follow the RPA frequency all the way to zero, due to numerical instabilities, although the trend 
is clear. Though we do not plot it, we also get a corresponding eigenvalue of the the stability matrix diving to zero. 

As the RPA frequency dives to zero, the corresponding magnitude of the hole-particle amplitude, |1a| 2 = 
Smi I^"™,a| 2 , increases dramatically. This we plot in the lower half of Fig. ©. As discussed in a previous sec- 
tion, it is this increase in \Y\ which causes the correlation energy to take on unphysically large values. Although we 
do not show it, we have also computed the expectation value for various operators, such as the Q ■ Q operator; while 
the HF contribution is stable, the RPA correlation correction Q also shows unphysically large values. 

We also considered 32 S, that is, 8 valence protons and 8 valence neutrons, and drove the 0c? 3 / 2 single-particle 
energy up and down. This case was very interesting because we could get, by adjusting A, spherical, prolate, and 
triaxial solutions. Fig. ([J} shows the shell-model (SM), Hartree-Fock, and HF+RPA energies for the ground state, 
while Fig. ([5]) shows the RPA frequencies and |1a| 2 - While a J— 2 (degeneracy = 5) eigenfrequency is falling, a J=4 
(degenercy = 9) mode actually falls below it, suggesting that ultimately it is a hexadecupole mode becomes unstable 
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FIG. 3: (Color online) For 28 Si. Upper panel: low-lying RPA frequencies Q for spherical (dashed) and deformed (solid) HF 
states. Lower panel: |Ya| 2 corresponding to the RPA frequencies in the upper panel. In both cases the degeneracy is given in 
parentheses; for the spherical case the degeneracy = 2 J + 1 where J is the angular momentum of the RPA mode. 
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FIG. 4: (Color online) Similar to Fig. ((2} but for 32 S; here A is the change in the 0^3/2 energy. 

first. The prolate state has both degenerate (time-reversed pairs) and non-degenerate states; the triaxial state, which 
has broken all possible symmetries, has no degenerate RPA modes. 

We reiterate: the key idea here is that the coexistence of locally stable phases leads to a first-order transition. The 
stable phases coexist easily because the parity-conserving mode-here quadrupole-means one can have cubic as well 
as quartic terms in the energy landscape. When we have an odd-parity mode, as illustrated in the next session, cubic 
terms are suppressed and one gets a second-order transition. 
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FIG. 5: (Color online) Similar to Fig. Q but for 32 S; here A is the change in the 0d 3/ / 2 energy. 

B. Example of a second-order transition: parity-mixing 

Second-order transitions lead to collapse of RPA. The classic example is the Lipkin model, which in its original 
form has a conserved parity (the Lipkin model is a two-level system, and the original Lipkin interaction could only 
promote or demote two particles at a time, thus providing an parity-like symmetry: either an even or an odd number 
of particles in the upper level); the transition of the HF state in the Lipkin model was from an exact parity state, 
with all particles in the lower level, to a mixed parity state. 

For our examples of second-order transitions, we consider a case involving shells of opposite parity. We fix the Os 
shell to be a closed core and have as the valence space 0p 3 / 2 , Opi/ 2 and 0d 5 / 2 , lsi/2- We look at 12 C, with 4 valence 
protons and 4 valence neutrons, without any truncations on the many-body space. (The only reason we leave out the 
0d 3 /2 orbit is to make full shell- model calculations tractable; SHERPA can easily handle the full p-sd space, but the 
HF+RPA results look very similar to what we present here. We also have looked at 16 and 20 Ne in similar spaces 
and get similar results.) We use the Cohen-Kurath (CK) matrix elements in the Op shell fl3l, the USD interaction [l2j 
in the 0d 5 / 2 -lsi/ 2 space, and the Millener-Kurath (MK) p-sd cross-shell matrix elements [141]. Within the p and sd 
spaces we use the original spacing of the single-particle energies for the CK and USD interactions, respectively, but 
then shift the sd single-particle energies up or down relative to the p-shell single particle energies by an amount A; 
we define A = where we get the first 3~ state at approximately 6.1 MeV above the ground state. The rest of the 
spectrum, in particular the first excited + state, is not very good, but the idea is to have a non-trivial model, not 
exact reproduction of the spectrum. 

(As a side note, this model space is not translationally invariant, and so we do not get zero-frequency modes from 
broken translational invariance.) 

Fig. [6] compares the exact SM ground state energy with HF and HF+RPA. Here the HF+RPA energy dives, or 
"collapses" at the transition point. A more detailed look is in Fig. [7| which plots the RPA frequencies. 

In Fig. [7] we also show the RPA eigenfrequencies and the magnitude of \Y\\. The exact-parity HF states are oblate, 
so that the RPA frequencies come in degenerate pairs from time-reversal symmetry. The mixed-parity HF states are 
triaxial, breaking time-reversal degeneracy. 

V. CONCLUSIONS 

We have applied Hartree-Fock plus the random phase approximation in the framework of the interacting shell 
model, with complicated, realistic forces. By changing the single-particle energies we could drive the Hartree-Fock 
solution between symmetry conserving and symmetry non-conserving states. In accordance with Thouless' original, 
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FIG. 6: (Color online) Ground state energy of 12 C in the 0pi/2-0pi/2-0d 5/ / 2 -lsi / / 2 space, calculated in the full interacting shell 
model (SM) (solid line), Hartree-Fock (HF) (dotted line), and Hartree-Fock plus random phase approximation (RPA) (dashed 
line). Here A is added to the lsi/2 and 0d 5 /2 single-particle energies; A = puts the negative parity states in the SM at 
approximately the correct location. We show where the HF states and corresponding HF+RPA states have either exact parity 
(and also spherical symmetry) or mixed parity (and break rotational invariance as well) 



8 



111,111, 




i i i i i , i i i 

(2) ~_ 


(1) 






(1) ~~ ~— - 

E 1 1 1 1 1 1 1- 


— i — i — Y—\ — i — \— 




- mixed parity 




H — — 1 — 1 — 1 — — 1 — 1 — \—, 
exact parity _ 


(1) — — — 




— ___ 


r 





(2) " 

1 



-6 -5 -4 -3 -2 -1 

A (MeV) 



FIG. 7: For 12 C. Upper panel: low-lying RPA frequencies Q for exact-parity and parity-mixing HF states. Lower panel: \Y\\ 2 
corresponding to the RPA frequencies in the upper panel. In both cases the degeneracy is given in parentheses; the exact- 
parity HF states are oblate, allowing for time-reversed pairs, but for parity-mixing the HF state becomes triaxial, breaking the 
time-reversal degeneracy. 
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and oft-forgotten, analysis, first-order transitions are associated with even-parity modes, such as the quadrupole mode, 
and do not display, at the transition point, the infamous collapse of RPA. Instead one only obtains the collapse of 
RPA in second-order transitions, associated with odd-parity modes. The latter are of course more serious, but it is 
useful to keep in mind that not all state or phase transitions automatically lead to the collapse of RPA. 

The U.S. Department of Energy supported this investigation through grants DE-FG02-96ER40985 and DE-FC02- 
07ER41457. 
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